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Resumen 


Ramos-Arzola, L. C., Cabrera-Estupiñán, E., Molina-Pérez, 
D., Hernández-Valdés, A. O., ££ Marón-Domínguez, D. E. 
(mayo-junio, 2017). Modelo para la optimización del costo 
de operación de un campo de pozos en acuíferos. Tecnología 
y Ciencias del Agua, 8(3), 39-53. 


El componente de mayor costo en la producción de agua 
subterránea en un sistema de suministro (una vez que los 
pozos están construidos) es el bombeo de agua desde los 
pozos. En la presente contribución se propone un nuevo 
módulo en MADA (modelo de administración de acuíferos), 
que se basa en la inclusión de un modelo de simulación- 
optimización para la minimización del costo de operación 
en acuíferos. La propuesta utiliza un modelo de simulación 
bidimensional del flujo del agua subterránea (AQUIMPE) 
basado en el método de los elementos finitos. Se emplea 
el enfoque matriz respuesta para vincular AQUIMPE con 
un problema de programación cuadrática, que se resuelve 
empleando la función quadprog del asistente matemático 
MATLAB. Se aplica el modelo en la administración de la 
explotación del acuífero Cuenca Sur y se obtiene una política 
de explotación que permite un 15% de ahorro del costo 
de bombeo en relación con el costo de la explotación real 
del año 2007. Además, el patrón de explotación obtenido 
satisface un conjunto de restricciones de demanda, de 
niveles en el acuífero y de capacidad instalada en cada pozo 
de explotación. 


Palabras clave: modelación del agua subterránea, modelo de 
simulación-optimización, enfoque matriz respuesta, MADA, 
AQUIMPE, programación cuadrática, Cuenca Sur de La 
Habana. 


Abstract 


Ramos-Arzola, L. C., Cabrera-Estupiñán, E., Molina-Pérez, D., 
Hernández-Valdés, A. O., € Marón-Domínguez, D. E. (May-June, 
2017). Model for optimization of the operating cost in aquifers. 
Water Technology and Sciences (in Spanish), 8(3), 39-53. 


The major cost component in the water production of groundwater 
supply system (once the wells are built) is pumping water from 
wells. In this contribution a new module in MADA (groundwater 
management model) based on the inclusion of a simulation- 
optimization model to minimize the operating cost in the aquifers 
exploitation is proposed. The proposal uses a two-dimensional 
simulation model of groundwater flow (AQUIMPE) based on the 
finite element method. The response matrix approach is used to 
link AQUIMPE with a quadratic programming problem, which is 
solved using the function quadprog of the mathematical assistant 
MATLAB. The model is applied in the administration of the 
exploitation in the Cuenca Sur aquifer, an exploitation policy is 
obtained that allows to save around a 15% ofpumping costs relative 
to the cost of the actual operation in the year 2007. In addition, the 
obtained exploitation pattern satisfies a set of constraints of water 
demand, aquifer levels and installed capacity in each operating well. 


Keywords: Groundwater modeling, simulation-optimization 
model, response matrix approach, MADA, AQUIMPE, quadratic 
programming, Cuenca Sur de La Habana. 
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Introducción 


Casi la mitad de la población mundial se abas- 
tece de agua potable por medio de fuentes sub- 
terráneas y se estima que 20% de los acuíferos 
está siendo sobreexplotado (WWAP, 2015). La 
enorme demanda sobre los recursos de agua 
subterránea (GWP, 2013), así como los efectos 
del cambio climático sobre este recurso motivan 
el desarrollo y empleo de técnicas que permitan 
administrar de forma racional el recurso hi- 
dráulico. Un problema bastante frecuente en la 
administración de acuíferos radica en obtener la 
política de explotación que garantiza el mínimo 
costo de bombeo de un campo de pozos en un 
sistema acuífero. 

Los modelos de simulación-optimización 
(MSO) o modelos de administración permiten 
encontrar alternativas para un uso sostenible del 
agua subterránea (Gorelick $: Zheng, 2015), re- 
solviendo un problema de optimización que se 
formula adecuadamente con objetivos definidos, 
y restricciones físicas y administrativas (Datta 
£z Kourakos, 2015). Para alcanzar una política 
de administración óptima, los MSO incorporan 
el modelo de simulación del agua subterránea 
como un conjunto de restricciones dentro del 
modelo de optimización (Yeh, 2015). El enfoque 
matriz respuesta (EMR) y el enfoque embebido 
(EE) son los dos métodos por lo general usados 
para incorporar el modelo de simulación dentro 


La energía requerida en la producción de 
agua subterránea es a menudo el mayor com- 
ponente energético y, por tanto, el mayor costo 
en un sistema de suministro de agua (una vez 
que los pozos están construidos), teniendo un 
mayor peso cuando se presentan bajos niveles 
del agua subterránea, producto de un bombeo 
sostenido (Ahlfeld $ Laverty, 2015). Para mi- 
nimizar el costo de operación o de bombeo es 
necesario optimizar una función que involucra 
el producto del caudal y la carga hidráulica, 
siendo esta última una función del caudal, por 
lo que es imposible abordar este problema con 
las técnicas de programación lineal y se requiere 
el uso de métodos de programación no lineal. 
En algunos casos, la carga hidráulica puede ser 
aproximada como una función lineal del caudal 
y, por tanto, el producto del caudal y la carga se 
convierte en una función cuadrática del caudal 
(Ahlfeld € Laverty, 2015). Ahlfeld y Laverty 
(2011, 2015) han obtenido una solución analí- 
tica mediante un MSO que garantiza un costo 
mínimo de bombeo. Sin embargo, esta solución 
descansa en las suposiciones de linealidad del 
sistema, condiciones de flujo estacionario y un 
valor adecuado de demanda externa. 

En Cuba se dispone de un modelo matemá- 
tico para la simulación del flujo subterráneo, 
AQUIMPE (Martínez, 1989), con amplias apli- 
caciones dentro y fuera del territorio nacional. 


Esta tecnología ha sido dotada con un modelo 


3 
pl 


del modelo de optimización (Gorelick, 1983). El 
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de administración de la explotación de acuíferos 


3, 
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primero usa el modelo de simulación de forma 
externa para obtener la respuesta de la carga 
hidráulica producto del bombeo y el segundo 
incorpora de forma explícita las ecuaciones del 
modelo de simulación dentro del modelo de op- 
timización (Ayvaz, 2009; Cabrera, 2009). Ambos 
enfoques han sido utilizados en la solución de 
varios problemas de administración, entre los 
que se encuentran maximización del rendimien- 
to de acuíferos, desarrollo de políticas para el 
uso conjunto de recursos de agua subterránea 
y agua superficial, protección de acuíferos cos- 
teros y minimización del costo de operación de 
un campo de pozos, entre otros (Gorelick, 1983; 
Ahlfeld $ Heidari, 1994; Wagner, 1995; Singh, 
2014; Gorelick $: Zheng, 2015; Yeh, 2015). 


llamado MADA (Cabrera € Dilla, 2011), que uti- 
liza el EE para acoplar AQUIMPE dentro de un 
problema de programación lineal. Por su propia 
concepción, MADA no es capaz de resolver la 
función no lineal que minimiza el costo de bom- 
beo; además, al utilizar el EE, considera los gas- 
tos y las cargas como variables de decisión del 
problema de optimización. Este último aspecto 
provoca un incremento significativo en la carga 
computacional en la medida que aumentan las 
dimensiones del problema (Taghavi, Howitt, dz 
Mariño, 1994; Yeh, 2015). 

En la presente contribución se crea un 
modelo matemático para resolver el problema 
de optimización que minimiza el costo de 
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operación de un campo de pozos destinado a 
captar agua subterránea. Dicha herramienta 
viene a conformar un nuevo módulo del modelo 
de administración del agua subterránea MADA. 
Para la concepción del nuevo modelo se utiliza 
el EMR para vincular AQUIMPE dentro de un 
problema de programación cuadrática que es 
resuelto mediante la función quadprog de MAT- 
LAB. Finalmente, el modelo se aplica al acuífero 
de la Cuenca Sur de La Habana, obteniéndose 
el patrón óptimo de explotación que garantiza 
el mínimo costo de bombeo, que representa un 
15% de ahorro en relación con el costo de la 
explotación real del año 2007. Además, el patrón 
de explotación obtenido satisface un conjunto 
de restricciones de demanda, de niveles en el 
acuífero y de capacidad instalada en cada pozo 
de explotación. Se demuestra que el plan de 
extracción del acuífero puede verse afectado 
por el efecto de varios años secos continuados. 


Materiales y métodos 
Modelo de simulación 


El modelo matemático AQUIMPE (Martínez, 
1989) permite simular el flujo impermanente 
bidimensional del agua subterránea a escala 
regional tanto en acuíferos libres como confi- 
nados, resolviendo numéricamente la siguiente 
ecuación: 


rn pH Ly 0 
OXxX 0x) oy dy ot 


sujeta a la condición inicial y las condiciones de 
borde que se muestran a continuación: 


h(x,y,0)=f.(x,y), (1,y) EQ (2) 
h(x,y,t)=f (2,4), (xy) €00, (3) 
TA) (uy) 500, (a 


donde h es la carga hidráulica; T, la transmisi- 

vidad del acuífero; S, el coeficiente de almace- 

namiento del acuífero; W, el término fuente o 

sumidero; t, el tiempo; x, y, las variables espacia- 

les; O, región de flujo; 9%, contorno de la región 
o) 


de flujo (92 = 4Q, U IQ,); rn derivada normal 
n 


al contorno, y f, f,, f, son funciones conocidas. 
AQUIMPE resuelve la ecuación (1) mediante 
el método de los elementos finitos con trián- 


gulo cuadrático y aplicando la aproximación 
de Galerkin (Cabrera éz Dilla, 2011). También 
permite simular la interacción del acuífero con 
un cuerpo de agua superficial (lago o embalse) 
siempre que exista comunicación entre ambos y 
posee un módulo para la calibración automática 
de las propiedades hidrogeológicas del acuífero 
(Gómez, Cabrera, € Garrido, 2009). 


Modelo de programación cuadrática 


El asistente matemático MATLAB está dotado 
con la función quadprog, que permite resolver 
un problema de programación cuadrática de 
acuerdo con la siguiente estructura: 


mín, 50H x+f” x 


em x= eg (5) 


Aye =b,, 


lb<x<ub 


donde x es el vector de las variables de decisión; 
H, la matriz hessiana de la función objetivo; f, un 
vector que representa los términos lineales de 
la función objetivo; A,,, " la matriz de las restric- 
ciones lineales de desigualdad; ¿0 el vector de 
los términos independientes de las restricciones 
lineales de desigualdad; A,,, la matriz de las res- 
tricciones lineales de igualdad; Do el vector de 
los términos independientes de las restricciones 
lineales de igualdad; los vectores lb y ub son las 
cotas inferiores y superiores, respectivamente, 
de las variables de decisión. 

La función quadprog es usada para resolver 
el modelo de administración propuesto y su 
sintaxis puede presentar la siguiente estructura: 
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[xopt, foal] 

=quadprog (Hf at dog! sd dor lb, ub) (6) 
donde xopt puede ser un punto de mínimo local 
para problemas no convexos y para problemas 
convexos ser un punto de mínimo global. El 
argumento de salida fvual es el valor mínimo 
local o global de la función objetivo. 

La matriz hessiana H (matriz de las segundas 
derivadas de la función objetivo) también puede 
ser utilizada como un medio para conocer si la 
solución encontrada es un óptimo global. Si H 
tiene la diagonal predominante y todos los ele- 
mentos de la diagonal principal son positivos, se 
puede decir que la función objetivo es convexa. 
Entonces, si las restricciones del problema de 
optimización son funciones lineales y existe la 
región factible de soluciones, todo el problema 
de programación cuadrática es convexo y el 
mínimo es global (Ahlfeld £ Mulligan, 2000; 
Ramos, Cabrera, Marón, € Hernández, 2016). 


Enfoque matriz respuesta 


En este enfoque, la ecuación (1), junto con las 
condiciones inicial y de borde, ecuaciones (2) 
y (3), son usadas para determinar la respuesta 
unitaria de la carga hidráulica, producto de 
un bombeo unitario en ubicaciones específicas 
dentro de la región de flujo (Yeh, 2015). En pro- 
blemas estacionarios, el EMR utiliza el principio 
de superposición en el espacio (Ramos et al., 
2016), y para problemas no estacionarios em- 
plea superposición en el tiempo y en el espacio 
(Psilovikos $: Tzimopoulos, 2004). El arreglo 
de todas las respuestas unitarias o coeficientes 
respuesta permite obtener la matriz respuesta 
del acuífero, la cual es introducida en el modelo 
de programación como un sustituto del modelo 
de simulación. 

Ahlfeld y Laverty (2011) señalan que cuando 
la ecuación de flujo del agua subterránea es 
lineal y las condiciones de borde no dependen 
de la carga hidráulica, la carga hidráulica es 
una función lineal del caudal y esta función se 
formula según la ecuación (7): 
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donde Q es el caudal; »%, la carga hidráulica 
cuando no se bombea en los puntos a optimizar 
el gasto; 9h /9Q, los coeficientes respuesta del 
acuífero; T, el número de periodos de admi- 
nistración; n, el número de puntos a optimizar 
el gasto; í, un punto donde se desea controlar 
la carga hidráulica, y ¡ es un punto donde se 
optimiza el gasto. 

El cálculo de los coeficientes respuesta se 
realiza mediante la ecuación (8), consistiendo 
en aproximar las derivadas parciales 9h /9Q por 
una fórmula en diferencia finita adelantada; este 
procedimiento también se denomina método de 
la perturbación (Ahlfeld, Barlow, £ Mulligan, 
2005): 


0h, , - h, ¿(Q,0) hi, 
Qs AQ; 


(8) 


donde AQ es una perturbación que se coloca 
en el punto ¡ a optimizar; Q, es un vector de 
gastos donde todas sus componentes son nulas, 
excepto en la componente j, donde es igual a 
AQ, y M(Q,,) es la carga en el punto donde se 
controla la carga luego de aplicada la perturba- 
ción. Es bueno señalar que como la respuesta de 
la carga al caudal es lineal, la derivada parcial 
9h /0Q es constante y la expresión (8) sólo está 
sujeta a errores de redondeo. Grava et al. (2015) 
subrayan el hecho de que para sistemas lineales, 
los coeficientes respuesta son independientes 
del valor de la perturbación. 

Para la obtención de todos los coeficientes 
respuesta y, por tanto, la obtención de la matriz 
respuesta, son necesarias n + 1 corridas del mo- 
delo de simulación (en este caso, AQUIMPE). En 
la primera corrida se considera que no existen 
extracciones en los puntos a optimizar, y se 
obtienen las cargas h” en todos los nodos y en 
todos los tiempos donde se controlará la carga. 
Luego, en las restantes n corridas se determinan 
las cargas h(Q m0) colocando la perturbación AQ 
en cada punto a optimizar. 
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La matriz respuesta, como se verá más ade- 
lante, permite transformar el problema de op- 
timización en un problema donde las variables 
de decisión son solamente los caudales en los 
nodos a optimizar el gasto. Este aspecto es muy 
conveniente cuando se administra un acuífero 
en régimen no estacionario, pues se reduce de 
forma considerable la dimensión del problema 
de optimización, en comparación con el EE, 
donde tanto las cargas como los caudales son 
variables de decisión. 


Formulación del modelo de administración 
Función objetivo 


El modelo MADA (Cabrera éz Dilla, 2011) con- 
sidera tres funciones objetivo lineales: (a) maxi- 
mizar las cargas, (b) minimizar los abatimientos 
y (c) maximizar el gasto de extracción. Además, 
MADA utiliza el EE y considera como variables 
de decisión los caudales de explotación en los 
nodos que se quieran optimizar en cada interva- 
lo de tiempo, Q,,, y también las cargas hidráu- 
licas en todos los nodos de la discretización del 
acuífero en cada intervalo de tiempo. En este 
trabajo se incorpora una cuarta función objetivo 
(9): minimizar el costo de bombeo que, como se 
dijo anteriormente, es una función cuadrática 
para el caso de sistemas lineales, y además se 
emplea el EMR para vincular AQUIMPE con 
un problema de programación cuadrática. El 
empleo del EMR permite considerar solamente 
como variables de decisión los caudales de 
explotación en los nodos a optimizar. 


T on 


minimizar F( 162 h) )-S CO; (H, -%,) (9) 


k=1 j=1 


donde H es la distancia entre la superficie del 
terreno y el plano de referencia en cada pozo 
(cota de la superficie del terreno), y c es un 
coeficiente que incluye términos para convertir 
de potencia a costo y toma la forma: 


Cc=cC,Aty (10) 


donde c, es el costo unitario de la potencia eléc- 
trica; At, la duración del bombeo, y y es el peso 
específico del agua. Según Ahlfeld y Laverty 
(2015), el coeficiente c también puede considerar 
la eficiencia de la bomba, factores de pérdidas 
por fricción y factores de conversión de unida- 
des. En la figura 1 se describen los términos de 
la ecuación (9). 

Mediante el EMR, la ecuación (9) puede 
ser escrita en términos de Q solamente como 
variables de decisión. Para ello se sustituye (7) 
en (9), realizando un cambio en los contadores 
de (7), y luego de sencillas manipulaciones se 
obtiene la siguiente expresión: 


minimizar F (Q) 


«ES Jo, [FI ico. Jo.) 


k=1 ¡=1 t=l m=1 


La función (11) debe ser derivada dos veces 
respecto a las variables de decisión Q para obte- 
ner la matriz hessiana necesaria para el modelo 
de programación cuadrática anteriormente pre- 
sentado. Luego de la derivación, los elementos 
de la matriz hessiana de F toman la forma: 


a [a ]_ Mn, (2) 
Or. Qu) “130, a 


Ahlfeld y Mulligan (2000) plantean que 
en la práctica los mayores valores de los ele- 


mentos de la matriz hessiana ocurren en la 
diagonal principal. Esto se debe a que el mayor 
abatimiento ocurre en el propio punto donde 
se bombea y, por tanto, la derivada respecto al 
gasto en esos puntos es grande (en términos 
absolutos, ya que 0h /0Q < 0). Entonces es bas- 
tante probable que la matriz hessiana tenga la 
diagonal predominante (y positiva) y F sea una 
función convexa, garantizando que el óptimo 
sea global. 


Restricciones 


Las restricciones de MADA pueden ser de tres 
tipos: (a) restricción de demanda, (b) restricción 


Tecnología y Ciencias del Agua, vol. VU, núm. 3, mayo-junio de 2017, pp. 39-53 


I55N 0187-8336 + 


Ramos-Arzola et al., Modelo para la optimización del costo de operación de un campo de pozos en acuíferos 


y Superficie del terreno Q 
/ Superficie piezométrica inicial 
sa erficie piezométrica 
Pp P 
H 
0 
Impermeable h h 
S ASA 
== 
7] 
E] 
E] 
rr 
LL] 
LL] 
Acuífero ] Datum 
= | 
e id, di a, ¡ao = a A AAA A 

E] 
A] 
E] 
= 
LS) 

Impermeable / 


Figura 1. Representación esquemática de los términos de la función que minimiza el costo de bombeo. 


de capacidad instalada, y (c) restricción de nive- En el caso de las restricciones de carga, 

les o carga hidráulica. Las dos primeras, ecuacio- ecuaciones (15) y (16), es necesario una trans- 

nes (13) y (14), respectivamente, son formuladas formación en términos de Q mediante la matriz 

en términos de los gastos que se optimizan y, respuesta. 

por tanto, su incorporación en el problema de 

programación cuadrática es sencilla. h.,=h ¿mín, k =1,..,T (15) 
II (13) h,, Ss h,¿máx, k=1,...,T (16) 
j-l 


Similar a la manera de transformar la fun- 
 _mínsQ,, <Q, máx, k=1...,T (14) ción objetivo (9), se sustituye la expresión (7) en 
Ñ : Ñ (15) y (16), y luego de acomodar se obtienen dos 
tricci función del iable d 
donde D, es la demanda prefijada en el periodo pe m ma O ÓN 
E o ñ ñ decisión Q, ecuaciones (17) y (18): 
de administración k, y Q, mín y Q,,máx son los 
caudales mínimo y máximo posible a extraer del 


SE 0h, 
pozo ¡en el periodo k. 22 Cn =h, mín-h0,,k=1,...,T (17) 
t=1 ¡=1 j. 
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T 
> =Q ;, Sh, ; máx- eL Tí (18) 
t=1 j=1 j,t 


Es conveniente señalar que las restricciones 
de carga no tienen por qué estar solamente aso- 
ciadas con puntos donde se desea optimizar el 
caudal, sino que pueden plantearse en cualquier 
lugar de la región de flujo. 


Modelo de administración 
El modelo de administración propuesto consiste 


en la formulación del problema de programa- 
ción cuadrática siguiente: 


minimizar 


E 


k=1 j=1 t=l m=1 


sujeta a: 


ik 0 
Qs: >h, , mín— ho , k=1l...,T 
e. * 


h 
3 O, sh,,máx-h?,, k=1,...,T 
an Ra * 


compuesto por las ecuaciones (11), (17), (18), 
(13) y (14), respectivamente. Este problema de 
programación cuadrática es resuelto mediante 
la función quadprog de MATLAB. El modelo pro- 
puesto sólo es válido para acuíferos confinados 
(o libres, que puedan ser tratados como confi- 
nados), flujo bidimensional, régimen imperma- 
nente y condiciones de borde no dependientes 
de la carga hidráulica. La secuencia de pasos 
para implementar el modelo de administración 
se muestra en la figura 2 y luego se aplica al 
acuífero de la Cuenca Sur de La Habana. 


Ya se ha expuesto que en la concepción 
inicial de MADA se utiliza el EE. Es conocido 
que con el EE se pueden presentar problemas 
de convergencia cuando se trabajan problemas 
de gran escala y complejidad (Taghavi et al., 
1994; Yeh, 2015). Sin embargo, en el EMR, los 
gastos son las únicas variables de decisión y las 
dimensiones del problema a resolver se reducen 
sustancialmente. 


Datos de entrada del 
modelo de simulación 
AQUIMPE 


Datos de entrada del 
modelo de administración 
MADA 


Obtención de la 
matriz respuesta 


(n + 1 corridas de AQUIMPE) 


Formulación del problema 
de programación cuadrática 


Solución del problema de 
programa cuadrática 
(función quadprog) 


Figura 2. Diagrama de flujo del modelo de administración 
propuesto. 
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Características de la zona de estudio provincia Artemisa, antigua provincia Habana 
(figura 3). Por sus características hidrogeo- 

El modelo de administración propuesto se lógicas, esta zona es de gran valor para la 
aplicó en el acuífero de la Cuenca Sur de La agricultura habanera y es también una de las 
Habana. El modelo existente de este acuífero principales fuentes de suministro de agua de 
se detalla en Cabrera y Dilla (2011) (en adelante la capital Ciudad de La Habana (Pérez, 2003). 
C£D). En ese trabajo se describen las principa- El medio geológico, principalmente compuesto 
les características hidrogeológicas del medio por calizas y dolomías del mioceno, posee una 
geológico, el modelo del acuífero y los datos gran heterogeneidad y anisotropía con dife- 
relacionados con los estímulos a los que está rentes grados de carsificación, además de una 
sometida esta fuente. Además, en CézD se apli- compleja composición química. Las principales 
ca MADA y se obtiene la distribución óptima hipótesis del modelo conceptual del acuífero son 
de caudales que maximiza el rendimiento del las siguientes: el acuífero puede ser simulado 
acuífero. A continuación se hace un resumen como confinado con un modelo bidimensional 
de los principales aspectos de este sistema no estacionario en planta, y el medio filtrante 
acuífero y de su modelo. puede simularse como si fuera un medio poroso 
La zona de estudio abarca un área de 830 km? isotrópico y heterogéneo, equivalente al medio 

y se encuentra enclavada en el sur de la real. En la figura 3 también se puede apreciar 


Zona de estudio 
a - 


Impermeable 


ro] 
200.000 400.000 600 000 metros 


Impermeable 
» 


Impermeable 


Carga constante, h = 0.2 m 


IA AP 
0 5 000 10 000 15 000 20 000 metros 


Figura 3. Ubicación de la zona de estudio, discretización del acuífero, condiciones de frontera y ubicación de los pozos de 
explotación del acueducto (círculos negros). 
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la discretización adoptada para la modelación 
del acuífero, en la cual se modificó la discreti- 
zación utilizada por CéfD (ver Cabrera, 2009). 
Esta nueva red de triángulos está compuesta por 
107 elementos o triángulos, y 234 nodos (entre 
principales y secundarios). 

En la figura 4a se observa la distribución 
de zonas hidrogeológicas obtenidas en una 
calibración previa de este acuífero (Llanusa, 
Martínez, Batista, € Hernández, 1993). Las 
conductividades hidráulicas se encuentran en 
un amplio rango desde 1 hasta 1 500 m/d, con 
espesores saturados entre 15 y 100 m. El valor 
del coeficiente de almacenamiento oscila de 0.1 
a 0.2. En el horizonte de planeamiento utilizado 
en el trabajo de CérD se consideró una recarga 
del acuífero proveniente de la lluvia del 75% 
de probabilidad, equivalente al comportamiento 
de un año medio seco, la cual fue distribuida 
espacial y volumétricamente en 12 grupos de 
triángulos (figura 4b). Para conocer a detalle la 
distribución de los valores de las propiedades 
hidrogeológicas y los valores de recarga del 
acuífero se puede consultar a CérD. 

En el contorno norte, el acuífero se encuentra 
en contacto con formaciones poco permeables, 
por lo que fue considerado como un contorno 
impermeable. Las fronteras este y oeste del 
modelo se asumieron como impermeables por 
hidrodinámica, al considerarse como líneas de 


corriente debido a que regionalmente existe un 
flujo del norte hacia el sur. Por último, el límite 
sur fue modelado como un contorno de descar- 
ga subterránea al mar, con carga conocida fija, 
con un valor de 0.2 m. En la figura 3 se pueden 
observar las condiciones de frontera descritas. 

Las principales extracciones realizadas al 
acuífero corresponden al acueducto (Cuenca 
Sur y provincia Habana) y a la agricultura. 
Según CéD, la suma de estas extracciones 
representa un 90% de la extracción global del 
acuífero. La extracción de agua por parte del 
acueducto de Cuenca Sur se realiza por medio 
de un campo de pozos compuesto por 19 pozos 
de explotación ubicados en la zona central del 
área estudiada (figura 3). Las características 
de interés de estos pozos del acueducto se 
muestran en el cuadro 1. 

En 2007, este campo de pozos extrajo un 
total de 97.58 hm”, que se distribuyeron (de 
forma espacial y temporal) según se muestra 
en el cuadro 2. En la última fila del cuadro 2 se 
puede apreciar la demanda mensual a la que fue 
sometido el acuífero por parte de este campo 
de pozos. 

Los pozos de extracción provenientes de la 
agricultura están dispersos por toda la zona de 
estudio y se dispone de muy poca información 
de las características de los mismos. Este aspec- 
to impidió que en este trabajo los pozos de la 


Figura 4. Zonas de propiedades hidrogeológicas (a), y zonas de igual comportamiento de la recarga (b). 
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Cuadro 1. Principales características de los pozos de explotación del acueducto Cuenca Sur. 


Cota de la superficie del terreno 


Nivel mínimo del agua subterránea 


Pozo de explotación Nodo la Ea) 

1 136 43.65 4.30 

2 146 41.15 4.36 

3 141 40.64 4.76 

4 142 35.48 3.60 

5 155 37.33 3.64 

6 84 32.61 3.04 

7 83 32.93 3.40 

8 72 32.07 3.39 

9 129 31.85 3.64 

10 111 31.39 3.52 

11 101 29.60 3.06 

12 124 10.29 1.91 

13 103 17.89 1.80 

14 89 18.23 1.90 

15 79 15.10 1.92 

16 65 13.07 2.01 

17 147 17.82 1.61 

18 158 17.03 1.80 

19 168 16.53 2.03 

Cuadro 2. Extracción real del campo de pozos de Cuenca Sur en el año 2007 (hm*/mes). 

Meses 
Pozo | Nodo | Enero | Febrero | Marzo | Abril | Mayo | Junio | Julio | Agosto | Septiembre | Octubre | Noviembre | Diciembre 

1 136 | 0.50 0.44 0.44 | 0.49 | 0.48 | 0.48 | 0.42 0.49 0.49 0.51 0.48 
2 146 0 0.28 0.44 | 0.47 | 0.48 | 0.51 | 0.44 0.49 0.49 0.51 0.48 
3 141 | 0.50 0.51 0.44 | 0.49 | 0.48 | 0.50 | 0.45 0.49 0.49 0.51 0.48 
o 4 142 | 0.50 0.51 0.44 | 0.49 | 0.48 | 0.51 | 0.46 0.48 0.48 0.51 0.48 
3 5 155 | 0.50 0.51 0.44 | 0.49 | 0.47 | 0.49 | 0.44 0.49 0.49 0.51 0.48 
=> 6 84 0.49 0.51 0.45 | 0.48 | 0.43 | 0.51 | 0.45 0.49 0.47 0.50 0.48 
E 7 83 0.49 0.27 0.45 | 0.48 | 0.44 | 0.50 | 0.44 0.49 0.47 0.50 0.48 
3 8 72 0.49 0.50 0.45 | 0.48 | 0.48 | 0.38 | 0.45 0.49 0.47 0.50 0.48 
E 9 129 | 0.50 0.50 0.45 | 0.48 | 0.49 | 0.51 | 0.46 0.49 0.46 0.45 0.40 
a] 10 | 111 | 0.50 0.50 0.45 | 0.48 | 0.49 | 0.51 | 0.47 0.49 0.43 0.40 0.47 
E 1 101 | 0.49 0.50 0.45 | 0.48 | 0.45 | 0.51 | 0.47 0.48 0.33 0,17 0.47 
E 12 | 124 | 0.49 0.39 0 0 0 0 0 0 0 0.18 0.50 0.48 
E 13 | 103 | 0.43 0.50 0.45 | 0.46 | 0.48 | 0.50 | 0.48 | 0.51 0.49 0.49 0.48 0.48 
5 14 89 0.50 0.06 0 0 0 0 0 0 0 0.04 0.50 0.45 
E 15 79 0.50 0.51 0.45 | 0.49 | 0.49 | 0.50 | 0.48 | 0.51 0.50 0.49 0.31 0.12 
= 16 65 0.50 0.23 0 0.35 | 0.49 | 0.49 | 0.48 | 0.51 0.50 0.50 0.50 0.24 
3 17 | 147 | 0.27 0.52 0.45 | 0.49 | 0.49 | 0.49 | 0.48 | 0.51 0.50 0.50 0.49 0.48 
E 18 | 158 | 0.21 0.38 0.45 | 0.41 | 0.48 | 0.49 | 0.46 | 0.51 0.50 0.49 0.49 0.48 
E 19 | 168 | 0.26 0.47 0.45 | 0.49 | 0.49 | 0.49 | 0.30 | 0.51 0.49 0.39 0.27 0.48 
8 Demanda | 8.12 8.09 7.15 | 8.00 | 8.09 | 8.37 | 7.63 | 8.63 8.35 8.15 8.61 8.39 
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agricultura fueran considerados como puntos 
a optimizar el caudal. Esto no quiere decir que 
las extracciones provenientes de los mismos no 
se consideraran, sino que se modelaron como 
puntos de extracción con caudal conocido, es 
decir, no entran en la optimización. 

La modificación a la discretización utilizada 
por CézD consistió en colocar un nodo de la 
malla en cada pozo de explotación que será op- 
timizado (pozos del acueducto de Cuenca Sur); 
de esta manera, se pueden reflejar de forma más 
adecuada los abatimientos que se generan en 
estos pozos, y tener un mejor control de la carga 
hidráulica en dichos puntos (ver la figura 3). 
Además, la adopción de esta nueva discretiza- 
ción sienta las bases para que en trabajos futuros 
se lleve una administración considerando los 
efectos locales que se producen en los pozos de 
bombeo (Cabrera €e Hernández, 2011). 

Otra modificación consistió en utilizar co- 
mo estado inicial de las cargas del acuífero el 
correspondiente con diciembre de 2006, donde 
se manifiestan valores de la carga de 0.2 a5m 
en la parte sur, y de 30 a 60 m en la parte norte. 

El Plan de Uso del Agua (PUA) propuesto 
por la Empresa de Aprovechamiento de los 
Recursos Hidráulicos de La Habana para el año 
1995 fue ajustado a un valor de 252.32 hm? /año. 
Este volumen, en 2007, se distribuyó en 97.58 
hm?, proveniente de la explotación del acueduc- 
to de Cuenca Sur; el resto, 154.74 hm”, pertenece 
a las extracciones de la agricultura y del acue- 
ducto de la provincia Habana. Por último, se 
consideró como demanda a suplir por el campo 
de pozos del acueducto de Cuenca Sur los va- 
lores que se muestran en la última fila (llamada 
demanda) del cuadro 2. 


Configuración para la aplicación de 
AQUIMPE y del modelo propuesto 


En el próximo apartado se realiza en AQUIMPE 
una primera simulación del acuífero sometién- 
dolo a la política de explotación real que se 
empleó en 2007. Siguiendo el criterio de ren- 
dimiento seguro de la escuela norteamericana, 


expuesto por CéD, se realiza otra simulación, 
repitiendo ese año hidrológico durante cinco 
años seguidos. Luego se aplica el modelo de ad- 
ministración propuesto para obtener el esquema 
de extracción que garantiza el mínimo costo de 
bombeo. Se consideran dos horizontes de pla- 
neamiento: el primero de un año y el segundo de 
cinco años, y se debe garantizar una demanda 
de 252.32 hm*/año. El valor del costo unitario 
de la potencia eléctrica se fijó en 0.3 $/kWh. 
Además, en la optimización se mantendrán 
constantes las extracciones provenientes de 
la agricultura y del acueducto de la provincia 
Habana. Con el esquema obtenido para el hori- 
zonte de planeamiento de un año se realiza una 
simulación durante cinco años seguidos. Todas 
las ejecuciones, tanto la simulación como la op- 
timización, consideran una recarga proveniente 
de la lluvia de 75% y las condiciones de frontera 
descritas. 


Resultados y discusión 


Aplicación de AQUIMPE con la explotación 
real realizada en el año 2007 


Se simuló el acuífero durante 2007 bajo las 
condiciones descritas. Esta simulación permitió 
obtener las cargas que ocurren en los pozos de 
explotación del acueducto de Cuenca Sur. Pos- 
teriormente se determinó el costo de bombeo en 
esos nodos, ecuación (9), arrojando un valor de 
$1.87-10% ese año. Repitiendo 2007 como patrón 
durante cinco años seguidos, en la figura 5 se 
puede observar la tendencia de regulación del 
acuífero. 

En la figura 5 también se puede apreciar 
que si se somete el acuífero a un periodo seco 
prolongado, los recursos explotables tienden al 
agotamiento, indicando que el plan de extrac- 
ción propuesto no garantiza la sustentabilidad 
de la fuente a largo plazo. Tampoco se puede 
asegurar que el patrón de explotación del cam- 
po de pozos del acueducto Cuenca Sur sea el 
de mínimo costo de bombeo para la demanda 
dada. 
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Figura 5. Gráfico de volumen almacenado en el acuífero durante cinco años para la explotación real del año 2007. 


Obtención y análisis del esquema de 
extracción de los pozos de Cuenca Sur que 
garantiza menor costo de operación en 2007 


Se considera un año como horizonte de pla- 
neamiento, con 12 periodos de administración 
(incremento mensual). Serán variables de deci- 
sión en el problema los caudales de los 19 pozos 
del acueducto de Cuenca Sur en los 12 meses 
del horizonte de planeamiento. El modelo de 
administración debe satisfacer las restricciones 
de niveles mínimos que se muestran en el 
cuadro 1; las restricciones de demanda que se 
muestran en el cuadro 2 (última fila) y el gasto 
en cada pozo estará limitado superiormente por 
0.532 hm*/mes. 

En esta primera ejecución se obtuvo la dis- 
tribución Óptima de caudales (cuadro 3) que 
garantiza el mínimo costo de bombeo y además 
satisface todas las restricciones anteriores. Con 
la solución obtenida, el costo de bombeo se 
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disminuyó hasta $1.59-10%, representando un 
ahorro de 15% ($280 000.00), en comparación 
con el costo de la explotación real del año 2007. 
En el cuadro 3 se puede observar cómo en mu- 
chos casos el modelo decide no extraer nada de 
algunos pozos de bombeo. Esto se debe a que 
dichos pozos están en zonas de baja transmisi- 
vidad y cuando el modelo evalúa un caudal en 
esos nodos, la carga desciende rápidamente y se 
viola la restricción (de carga mínima). 

Otro aspecto de interés de la solución obteni- 
da está en la naturaleza del óptimo encontrado. 
Esto se pudo comprobar, ya que la matriz he- 
ssiana del problema tiene la diagonal predomi- 
nante y positiva, entonces la función objetivo 
es convexa, y como se trata de un problema 
con restricciones lineales de desigualdad, todo 
el problema es convexo y el óptimo encontrado 
es global. 

Sin embargo, con la solución óptima obte- 
nida tampoco se garantiza que a largo plazo 
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Cuadro 3. Extracción óptima del campo de pozos de Cuenca Sur (hm? /mes). 


Meses 
Pozo Febrero Mayo Diciembre 
1 0 
2 0 
3 0 
4 0.53 
5 0.41 
6 0.53 
7 0.53 
8 0.53 
9 0.53 
10 0.53 
11 0.53 
12 0.53 
13 0.53 
14 0.53 
15 0.53 
16 0.53 
17 0.53 
18 0.53 | 0.53 0.53 
19 168 | 0.53 0.53 0.53 | 0.53 | 0.53 | 0.53 | 0.53 | 0.53 0.53 0.53 0.53 0.53 


(cinco años) el acuífero sea capaz de satisfacer 
la demanda sin violar las restricciones de nivel. 
Esto se puede evidenciar al observar la curva 
con triángulos en la figura 6, donde a partir 
del cuarto año los niveles en el pozo 5 (nodo 
155) comienzan a estar por debajo del nivel 
mínimo permisible; aunque, como se aprecia, 
con la explotación real los niveles descienden 
más rápidamente y se violan antes los límites 
de carga mínima. 


Obtención y análisis del esquema de 
extracción de los pozos de Cuenca Sur que 
garantiza menor costo de operación; horizonte 
de planeamiento: cinco años 


Para esta ejecución, el problema de optimi- 
zación constó de 1 140 variables de decisión 
(los gastos de los 19 pozos durante 60 meses) 
y 1 200 restricciones de desigualdad. La curva 
con cuadrados de la figura 6 evidencia que no 
fue posible obtener una solución factible y a 


partir del cuarto año el modelo fue incapaz 
de garantizar que los niveles en el acuífero se 
encontraran por encima de los niveles míni- 
mos. Esta situación se venía manifestando en 
las ejecuciones anteriores, donde se apreciaba 
que la demanda real a la que estaba sometida el 
acuífero no garantizaba la sustentabilidad de la 
fuente para un periodo seco prolongado. Para 
esta solución no factible se obtiene un costo de 
bombeo de $8.64-10*. 


Conclusiones 


Se crea un modelo matemático que permite mi- 
nimizar el costo de bombeo de la explotación de 
acuíferos. Esta herramienta viene a conformar 
un nuevo módulo del modelo de administración 
de acuíferos MADA. El nuevo modelo utiliza 
el enfoque matriz respuesta para vincular el 
modelo de simulación de acuíferos AQUIMPE 
dentro un problema de programación cuadrática 
que es resuelto mediante la función quadprog de 
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Figura 6. Carga mínima y cargas simuladas en el pozo 5 durante un periodo de cinco años. 


MATLAB. Además, se tienen en cuenta restric- 
ciones de demanda, nivel y capacidad instalada. 
El modelo propuesto sólo es válido en acuíferos 
hidráulicamente confinados. 

Se obtiene la distribución Óptima de la ex- 
plotación del campo de pozos del acueducto 
Cuenca Sur que garantiza el mínimo costo de 
bombeo, obteniéndose un ahorro de 15% en 
relación con el costo de la explotación real del 
año 2007 a la que fue sometida este acuífero. 
Además, la solución obtenida satisface todas las 
restricciones administrativas y físicas impuestas. 
De esta manera, queda demostrada la utilidad 
práctica del modelo propuesto. 

No fue posible obtener una solución factible 
en el caso de que el acuífero esté sometido a un 
periodo seco prolongado, pues el modelo fue 
incapaz de garantizar que los niveles en el acuí- 
fero se encontraran por encima de los niveles 
mínimos. Esto indica que el plan de extracción 
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del acuífero puede verse afectado por el efecto 
de varios años secos continuos. 
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